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We show that quasiparticles in a magnetic field of arbitrary strength H can be described by field 
dependent parameters. We illustrate this approach in the case of an Anderson impurity model and 
use the numerical renormalization group (NRG) to calculate the renormalized parameters for the 
levels with spin a, i d ,cr(H), resonance width A(H) and the effective local quasiparticle interaction 
U(H). In the Kondo or strong correlation limit of the model the progressive de-renormalization 
of the quasiparticles can be followed as the magnetic field is increased. The low temperature be- 
haviour, including the conductivity, in arbitrary magnetic field can be calculated in terms of the 
field dependent parameters using the renormalized perturbation expansion. Using the NRG the field 
dependence of the spectral density on higher scales is also calculated. 



INTRODUCTION 



II. QUASIPARTICLES AT T = WITH ZERO 
MAGNETIC FIELD 



Electrons in strongly correlated systems are particu- 
larly sensitive to the application of magnetic fields. One 
reason is that strong correlations are usually a conse- 
quence of the interaction of electrons with enhanced spin 
fluctuations, and these fluctuations couple strongly to a 
magnetic field. Another reason is that there is a low tem- 
perature scale T* (T* <C Tp) induced which plays the 
role of an effective Fermi temperature Tf. The effects 
of a magnetic field H in general depend on the ratio of 
the two energy scales /j-bH and k^T-p. In a weakly corre- 
lated metal (J-bH / 'k^Tp <C 1, but in a strongly correlated 
system the relevant ratio is usH/kBT*, which can be of 
order unity. This sensitivity means that a magnetic field 
is an important tool in the experimental investigation of 
strongly correlated metallic systems, such as magnetic 
impurities, quantum dots, heavy fermions and transition 
metal oxides. 



In this paper we concentrate mainly on the effects of 
a magnetic field on the quasiparticles in the Fermi liquid 
regime for various models of strongly correlated metals. 
In particular we develop an approach, based on a combi- 
nation of the numerical renormalization group (NRG) 1 ' 2 
and renormalized perturbation theory (RPT) 3,4 . This 
approach gives a comprehensive picture of the behaviour 
of quasiparticles in magnetic fields of arbitrary strength, 
such that we can follow the renormalization or de- 
renormalization of the quasiparticles as the magnetic field 
strength is changed. In this paper we concentrate on de- 
veloping the approach for the single impurity Anderson 
model (SIAM), as a model of magnetic impurities and 
quantum dots. This prepares the way for generalizing 
the approach to lattice models for transition metal ox- 
ides and heavy fermion materials, using dynamical mean 
field theory (DMFT)£, which will be the subject of a 
subsequent paper. 



The generic model for strongly correlated local sys- 
tems, such as magnetic impurities in a host metal or a 
quantum dot coupled to an electron reservoir, is the An- 
derson model 6 . The Hamiltonian for this model is 

Ham = 22 e d ^d\d a + Un d ^n d ,i (1) 
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where e dM = £d — ^gfJ-sH /2 is the energy of the localized 
level at an impurity site or quantum dot in a magnetic 
field H, U the interaction at this local site, and 14 the 
hybridization matrix element to a band of conduction 
electrons with energy e k - When U = the local level 
broadens into a resonance, corresponding to a localized 
quasi-bound state, whose width depends on the quantity 
A(lu) — irJ2k \Vk\ 2 o~{^ — £k)- It is usual to consider the 
case of a wide conduction band with a flat density of 
states where A(w) becomes independent of u> and can be 
taken as a constant A. 

The low energy behaviour of this model can be ex- 
pressed in terms of the renormalized quasiparticles of a 
local Fermi liquid, which is described by a renormalized 
version of the same model: 

Ham = Sddjd^ + U : n d ^n dd : 

a 

+ Y,Wl C k,a + V k *4 !<T d a ) + ]T Sk,J Ka C k ^ (2) 
k,a k,a 

where the colon brackets indicate that the expression 
within them must be normal-ordered. This Hamiltonian 
corresponds to the low energy fixed point of the Wilson 
numerical renormalization group transformation of the 
discrctized Anderson and Kondo models, with the lead- 
ing irrelevant termsiSi The advantage of describing 
the fixed point in this way, as a renormalized Anderson 
model rather than as a strong coupling fixed point of the 
Kondo model, even in the strong correlation or Kondo 
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limit, is that it clearly brings out the 1-1 correspondence 
of the low-lying single particle excitations with those of 
the non-interacting model£*S*^. Furthermore, it is ap- 
plicable in all parameter regimes, from weak to strong 
coupling and for all occupation values for the local site. 
The parameters id, A, and U define the quasiparticles of 
this renormalized model, and a simple direct procedure 
for calculating these parameters using the NRG has been 
given earliest. In terms of these parameters the Friedel 
sum rul o 1 ^ 11 , which gives the total occupation of the d- 
orbital at the impurity site rid for T = in the wide band 
limit, is 



1 2 f -i ( id 
rid = 1 tan — 

7T V A 



(3) 



In the wide band limit the renormalized parameters, 
id, A, and U, can be expressed as functions of two 
independent 'bare' parameters, which can be taken to 
be U/irA and Ed/VA. Some typical plots of the renor- 
malized parameters as a functions of these variables 
were presented in the earlier work£. An alternative way 
of presenting the results is in terms of just one of these 
variables, U/ttA, and the occupation of the impurity 
levels rid, which gives a global picture over the various 
regimes of the model. 

In figures 1(a), 1(b) and 1(c) we present results for 
the ratios, A/A, id/id and U/U, where e~d = £d + U/2, 
which give a measure of the degree of renormalization 
of these parameters for the 'bare' values of U/irA = 
0.5,1,1.5,2,2.5,3. In the empty (n d -> 0) and full 
(rid ~~ * 2) regimes all the parameters approach their un- 
renormalized values so these ratios tend to unity. The 
renormalizations are the most pronounced in the region 
rid — 1. In the almost localized (Kondo) regime, rid — 1, 
U/irA > 2, there is only one energy scale, the Kondo 
temperature 7k, such that £d — and U = 7rA = 47k. 

The impurity density of states for the non-interacting 
quasiparticles p~d(w) at T = is given by 



A/tt 



(w - £d) 2 + A 2 



(4) 



and in the localized regime corresponds to a Kondo reso- 
nance of half-width A = 47k/7r at the Fermi level, where 
the Kondo temperature is defined in terms of the impu- 
rity spin susceptibility %s at T = via 7k = (5Mb) 2 / '4Xs- 
As in the general Fermi liquid theory, the linear coeffi- 
cient of the specific heat for the impurity can be expressed 
simply in terms of the quasiparticle density of states at 
the Fermi level, and is given by 



7d 



2tt 2 , . 

-g-/5d(0). 



(5) 



Exact expressions for impurity spin and charge suscepti- 
bilities, Xs and Xci m terms of the renormalized param- 
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FIG. 1: Plots of the renormalized parameters for the 
Anderson model, (a) A/A, (b) idji& and (c) U/U, 
as a function of the occupation of the impurity site 
nd, for values of U/ttA = 0.5,1.0,1.5,2.0,2.5,3.0. 



ctcrs for T = are 
1 

Xs 



1 



2 Pd(0)(l + Up d (0)), xc - 2/5d(0)(l - Up d (0)). 

■:<ii 

These results can be derived using the renormalized per- 
turbation approach^ working simply only up to first or- 
der in U (see also the Appendix here). 

It is interesting to note that these Fermi liquid results 
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apply even in the limit Sd — > — oo, U — > oo, when the 
charge fluctuations are completely suppressed and the 
Kondo scale Tk — > 0. In this limit Upd(0) — > 1 and 
— ► <5(w) corresponding to the low energy excita- 
tions of an isolated spin. If the finite temperature Fermi 
distribution function is included in the calculations lead- 
ing to Xc and Xs in equation JSJ), then in this limit Xc - > 0, 
and Xs — > (5MB) 2 /4fcB7 1 , the Curie law for a localized 
magnetic moment. 

So far we have assumed a repulsive interaction U > 0, 
but the formulae are equally applicable to the case of an 
attractive interaction U < 0. In this case, however, ex- 
cept at the particle-hole symmetric point, the renormal- 
ized parameters behave quite differently as a function of 
\U\/tvA and rid. In figures 2(a), 2(b) and 2(c), we present 
results for the ratios, A/A, £~d/£~d, and U /U, for values 
of U/ttA = -(0.5, 1, 1.5, 2, 2.5, 3). 

In comparing the results in figure 2(a) (b) and (c) with 
the corresponding results in figure 1 (a)(b) and (c), it 
can be seen that, irrespective of the sign of U all the pa- 
rameters converge to their bare values in limits rid — ► 
and rid — > 2 and the values of A/ A and U/U are the 
same at the symmetric point rid = 1, they are remark- 
ably different elsewhere. In the range 0.5 < rid < 1.5, 
A/ A and U/U for U < are remarkably flatter than the 
corresponding results for U > 0, especially for large val- 
ues of \U\. The ratios £d/£~d and U /U for U < do not 
increase monotonically with increase of |rid — 1|, as they 
do in the corresponding U > results. The ratio U /U 
develops a sharp peak with values such that U > U and 
is completely different from the positive U counterpart. 
The ratio £d/£d is also completely different and is greater 
than, or equal to, unity over the whole range, developing 
a similar peak to that for U/U. At the particle-hole sym- 
metric point, rid = 1, £d = 0, and £~d = 0, so the ratio 
£d/£d is not defined. In the positive U case, the ratio 
tends to for large \U\ and to 2 in the negative U case. 

The main applications of the Anderson model to mag- 
netic impurities and quantum dots are for the positive 
U case, though the negative U model has some limited 
application as an effective model for some locally coupled 
electron-phonon systems, where a local attraction is in- 
duced through the exchange of a virtual phonor. 12 . How- 
ever, in the next section we shall exploit the negative U 
model here by using spin-isospin symmetry to transform 
it into a model with positive U with particle-hole sym- 
metry in the presence of a magnetic field. The restriction 
to the rid = 1 particle-hole symmetric case will not be a 
serious limitation because in the strong correlation/local 
moment limit the renormalized level rapidly approaches 
the Fermi level for large U, £d — > 0, and hence via the 
Friedel sum rule rid — > 1. In this particle- hole symmet- 
ric case there is only one energy scale for a given U/tvA, 
which we denote by T* and define by T* — n A/ 4; it is 
such that T* — > T K in the Kondo regime U/irA > 2. 
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FIG. 2: Plots of the renormalized parameters for the 
Anderson model, (a) A/A, (b) £~d/£d and (c) U/U, 
as a function of the occupation of the impurity site 
n d , for values of U/ttA = -(0.5,1.0,1.5,2.0,2.5,3.0). 



III. QUASIPARTICLES AT T = WITH 
ARBITRARY MAGNETIC FIELD H 

The spin-isospin transformation which formally elimi- 
nates the magnetic field term in the symmetric model is 
defined by 
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and 



„t 

-fc-p 



--kl 



(8) 



with the requirement e_fc = — with V^ fe = — Vfc. 
With this transformation for the symmetric model with 
[7 > and £d = in the presence of a magnetic field 
H at the impurity site gets mapped into the model 
with an interaction —U, and no magnetic field but with 
£d = = h. The spin/charge susceptibilities 

for the positive U symmetric model in a magnetic field 
are then given by the charge/spin susceptibilities of the 
model with a negative U. The magnetic field depen- 
dent parameters A(h), £d,a(h) and U(h), for the U > 
model then correspond to A(h), — cre&(h) and —U(h) for 
the corresponding negative U model with Ed = h. The 
Friedel sum rule is still applicable to each spin component 
from which the induced impurity magnetization M(h) at 
T = can be deduced, 



m(h) 



M(h) 
5MB 



2 (nd < T 



1 _! /£ d (/l)>| 

rid I ) = -tan -= , 

(9) 

from the two parameters £d(M and A(h) that character- 
ize the non-interacting quasiparticles. The quasiparticle 
density of states for the spin up and spin down electrons 
is given by 



Pd >a (uj,h) 



A(h)/w 



{u-ai d {h)) 2 + A 2 (h)- 



(10) 



and the field dependent spin and charge susceptibilities 
at T — are given by 

Xs(h) = ^p d {0,h)(l + U(h)p A {0,h)), 



of the effective interaction U(h), as the magnetization is 
reduced from the saturated value, m sat = 1/2. As the 
applied magnetic field is reduced from the regime h > U, 
spin fluctuations increase and enhance the effective inter- 
action U, as in the random phase approximation (RPA), 
above the bare value U (we will see this more explicitly 
later). As the magnetic field is further reduced the many- 
body correlations are increasingly effective in screening 
the impurity so that U(h) decreases from an enhanced 
value greater than U to a value 4Tk as h — ► when 
U > 2ttA. 




FIG. 3: The magnetic field dependence of the renormal- 
ized parameters A(h)/A, ed(h)/e~d (= fj(h)) and U(h)/U, 
and the Wilson ratio R(h), for the symmetric Anderson 
model with U/irA — 3.0 plotted on a logarithmic scale. 



Xc(h) = 2^(0, h)(l - U(h)p d (0, h)). (11) 

As /5d,cr(0, h) is independent of the spin state we can drop 
the spin index a. 

As id,<r(h) is entirely magnetic field driven it is con- 
venient to write it as id,a(h) = hfj(h), then 2hfj(h) is 
the Zeeman splitting of the impurity levels for the non- 
interacting quasiparticles. 

Using these results we can reinterpret the renormal- 
ized parameters in figures 2(a), (b) and (c) for the U < 
model as A(h)/A, fj(h) and U(h)/U for the U > model 
as functions of the variable 2m(h) + 1, which replaces the 
rid- The fact that £d/£d in figure 2(b) approaches the 
value 2 as rid ^ 1 in the strong correlation limit can now 
be interpreted as ij(h) — > 2 for h — > 0, which is equivalent 
to the Wilson ratio R = 47rx s /3(.g/JB) 2 7imp = 2 for the 
strongly renormalized quasiparticles in zero field, which 
is enhanced compared with the free electron value R = 1. 
In figure 2(c) the sharp rise in U/U on reducing rid from 
the value rid = 2, can be interpreted as the enhancement 



In figure 3 we give a more conventional plot of the 
renormalized parameters as a function of the natural log- 
arithm of the magnetic field, ln(/i/T*), for the strong 
coupling case U/nA = 3, where T* = ttA(0)/4 = 2.00 x 
10~ 2 = Tk, in agreement to within 0.5% with formula^ 
for T K for this model, T K = ^([/A/2)e- 7rC/ / 8A+TrA / 2l/ . 
We can follow the progressive de-renormalization of the 
quasiparticles as the strong correlation effects are sup- 
pressed as the magnetic field is increased. Initially the 
quasiparticle interaction U (h) increases and can reach 
values such that U(h) is greater than the bare interac- 
tion U. This does not imply, however, that the inter- 
action effects are becoming stronger. The effects of the 
interaction on the low energy scale depend upon the com- 
bination, U(h)pd(0, h), and pd(0, h) falls off rapidly with 
h as £d(h) moves away from the Fermi level. The Wilson 
ratio R(h) = 1 + U(h)pd(0, h) is a measure of this com- 
bination of factors and for the Kondo model it is known 
from Bethe ansatz calculations^ that R(h) — 2 is inde- 
pendent of h. This can be seen to be the case in the 
results for R(h) shown in figure 3 when the parameters 
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FIG. 4: The impurity magnetization m(h) for the sym- 
metric model with U/irA = 3.0, together with i?(ft)/4, 
where R(h) is the Wilson ratio, plotted as a function 
of the logarithm of the magnetic field. Also shown for 
comparison are the corresponding Bethe ansatz results^* 
for the field induced magnetization for the Kondo model. 



correspond to the localized or Kondo regime. The local- 
ized model, however, is only valid when the charge fluc- 
tuations are completely suppressed. For very large field 
values h > U local charge fluctuations can be induced by 
the magnetic field and, as this regime is approached, R(h) 
makes a crossover to the value R — 1 for non-interacting 
electrons. The combination U(h)p d (0,h) can be seen to 
decrease monotonically with increase of ft. 



Equation for the susceptibility Xs (ft) has a term in 
U(h). However, the susceptibility can also be derived by 
differentiating the expression @ for the magnetization 
which depends explicitly only on the variables e d (h) and 
A(ft). Hence the value of U(h) is not independent of 
the other two parameters and we can derive a relation 
between them, 



1 + U(h)p d (0,h) 



de d (h) e d (h) dA(h) 



dh 



A(h) dh 



(12) 



The proof that equation for the susceptibility is ex- 
act depends on a Ward identity, so the relation (|12[1 we 
have derived must be an alternative statement of this 
identity. In terms of fj(h) — e d (h)/h it becomes 



l + U(h)p d (0,h) =fj(h) + h 



dfj(h) hrjjh) dA(h) 
dh A(h) dh 



(13) 



For ft = it implies that fj(0) — i?(0), which can be 
seen in the results in figure 3. 

Mean field theory can also be interpreted in terms of 
renormalized parameters with e d (h) = ft + UrriMF(h), 



A(ft) = A, where mMF(ft) is the mean field magnetiza- 
tion. These parameters are substituted into equation JjJJ , 
and mMF(ft) is determined self-consistently. From equa- 
tion (|12|) we can deduce the corresponding value of U(h) 
which gives 



U(h) = 



U 



1 - Up d MF{Q,h)'' 



(14) 



where PdMF(w, ft) is the quasiparticle density of states 
with the mean field parameters. This result corresponds 
to the enhancement of the susceptibility that one finds 
from the random phase approximation. If the magnetic 
field is reduced from a large value ft > U then pdMF(0, ft) 
decreases and so U(h) increases. This is precisely what 
is seen in the large ft regime in the results in figure 3, as 
well as those shown earlier in figure 2 (b). 

In figure 4 we plot the magnetization derived from 
these parameters for U/ttA = 3 as a function of ln(ft/T*) 
using equation We also give the magnetization de- 
duced from the Bethe ansatz results for the Kondo model 
for comparison together with the value of i?(ft)/4. There 
is complete agreement with the results of the Kondo 
model, up to and just beyond the point at which local 
charge fluctuations are induced by the magnetic field, 
where R(h)/4 begins to decrease significantly from its 
strong correlation value 0.5. The ratio s d (h) / AQi) must 
be a universal function of ft/Tx in the Kondo regime. 
For ft > Xk we have the asymptotic form 1/2 — m(ft) ~ 
l/41n(ft/TK) 15 , so from equation ©, in this regime the 
ratio £d(ft)/A(ft) cx ln(ft/Tk). In figure 5 we plot the ra- 
tio £d(ft)/A(ft) against ln(ft/T*), for U/nA = 5, which is 
well in the Kondo regime so T* — Tk- It can be seen this 
ratio is proportional to ln(ft/Tk) for a significant range of 
magnetic field values for ft > Tk, before the effects of the 
field induced charge fluctuations take over. Due to the 
charge fluctuations, the approach to saturation is much 
more rapid for the Anderson model than for the Kondo 
model, once ft exceeds U, as can be seen clearly in figure 
4. 

A similar plot of the magnetization m(ft) against field 
on a logarithmic scale in shown in figure 6 for a weak 
correlation case U/ttA = 0.25. In this regime it is more 
appropriate to compare the results with mean field theory 
The mean field results for the same parameters are also 
given in figure 6 (dashed line) and can be seen to be in 
good agreement. Also shown is the Wilson ratio i?(ft)/4, 
which is only weakly enhanced at ft = 0, R(0) — 1.244, 
and makes a slow crossover to that for non-interacting 
electrons for ft ~ A. 

In figure 7 we give the magnetization m(ft) against 
ft/T* for the more realistically realizable magnetic field 
regime < ft < 2.5T* for a strong correlation case 
U /ttA — 3, using the renormalized parameters in equa- 
tion Points (x) corresponding to Bethe ansatz 
resultsi- for the Kondo model, are included for compar- 
ison and can be seen to be in complete agreement with 
those for the Anderson model in this regime. The re- 
sults for m(ft) for the weak coupling case U/irA — 0.25 
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FIG. 5: The ratio e d (h)/A(h) for U/nA = 5.0 plot- 
ted as a function of the logarithm of the magnetic field. 




Ln(h/T*) 



FIG. 6: The impurity magnetization m(h) for the 
symmetric model for a weak correlation case with 
U/nA = 0.25, together with R(h)/A, where R(h) is 
the Wilson ratio, plotted as a function of the loga- 
rithm of the magnetic field. Also shown for com- 
parison are the corresponding results (dashed line) for 
the magnetization calculated using mean field theory. 



are also shown in figure 7 together with the correspond- 
ing mean field theory results (dashed line). The T*s for 
the strong and weak correlation cases are very different, 
T*(3)/T*(0.25) = 0.0827, so the energy scales are very 
different but, relative to these scales, m(h) increases ini- 
tially more slowly with h in the weak correlation case 
but approaches saturation more rapidly. The mean field 
results U/nA = 0.25 also give a good approximation to 
m(h) for U/nA = 0.25 over this range. 

Also shown in figure 7 are the values of R(h)/A for the 
strong correlation case U/ttA = 3. It can be seen that 
there is a 5% reduction in R(h) over this range to h — 
2.5Tk, indicating some evidence of charge fluctuations 
beginning to contribute to the specific heat coefficient 



FIG. 7: The impurity magnetization m(h) and R(h)/4, 
where R(h) is the Wilson ratio, for the symmet- 
ric model with U/nA = 3.0 plotted as a function 
of the magnetic field h/T*, and compared with the 
Bethe ansatz results^ (crosses) for the Kondo model. 
Also shown are the results for m(h) calculated for 
the weak coupling case U/nA = 0.25 together with 
the corresponding mean field results (dashed line). 



7d(/i) for intermediate field values, but having little effect 
on the magnetization. The localized model gives R(h) = 
2 for all h, implies that U(h) = l/p^(0,h). From this 
result, and equations (jTU|> and ©, the ratio U(h)/nA(h) 
for the localized model can be expressed entirely in terms 
of the magnetization and is such that U(h)/wA(h) = 
1/ cos 2 (nm(h)). For h = 0, this corresponds to the strong 
correlation results U(0)/nA(0) = 1, as m(0) = 0, and 
for very large fields where m(h) — > 1/2 as h — > oo, it 
gives U(h)/nA(h) — > oo, corresponding to the fact that 
charge fluctuations can only be completely suppressed if 
U is infinite. In figure 8 we plot the ratio U(h)/nA(h) 
as a function of \n(h/T*) for U/nA = 3, 5, and compare 
the results with 1/ cos 2 (7rm(/i)), where m(h) is calculated 
from the Bethe ansatz results for the Kondo model. It 
can be seen that the large peak in U(h)/nA(h) develops 
as h increases, because it closely follows the result for the 
localized model 1/ cos 2 (7rm(/i)), and then falls back in 
very large fields to the unrenormalized ratio U/nA. The 
results for U/nA — 5 follow the form for the localized 
model 1/ cos 2 (7T77i(ft,)), implying R(h) = 2, for a greater 
range of h than for the case U/nA = 3, resulting in a 
more pronounced peak. 



IV. LOW TEMPERATURE 
THERMODYNAMICS IN AN ARBITRARY 
MAGNETIC FIELD 

We can generalize some of the results in the preceding 
section to calculate the leading temperature dependence 
in the presence of an arbitrary magnetic field. For in- 



7 




Ln(h/T») h/T* 



FIG. 8: The ratio U(h)/nA(h) is plotted as 
a function of \n(h/T*) for the cases U/irA — 
3, 5. The two results are compared with the re- 
sult 1/ cos 2 (Trm(h)) , where m(h) is calculated from 
the Bethe ansatz results^ for the localized model 
where all impurity charge fluctuations are suppressed. 



FIG. 9: The ratio — p' d (h) / p^(0) , where the prime 
indicates a derivative with respect to h/T*, is 
shown for U/irA = 3.0,0.5,0.0 as a function of 
h/T*. The dotted line shows the asymptotic re- 
sult as h — + 0, 3/171-^3/27^, for the Kondo model. 



stance, for the temperature dependence of the suscepti- 
bility and magnetization we can use the thermodynamic 
relation. 



d 2 C(T,H) ^ T f X (T,H) 



dH 2 



8T 2 



(15) 



where C(T, H) is the specific heat. Applying this to the 
impurity contribution and taking the limit T — > we find 



d 2 Xs(T,H) 



dT 2 



T=0 



d 2 7d (g) 
dH 2 ■ 



(16) 



We can use the result in equation (J5J to deduce the results 
for the T 2 dependence for Xs (T, h) , 



Xs(P,h)- X ,(T,h) 



7T 2 T 2 d 2 ~ Pd {Q,h) 
12 Xs (0,/i) dh 2 

2 



cJh) 



On integrating these results with respect to h we can 
derive a similar relation for the induced magnetization, 
M(T,h) = m(T,h)/(gp B ), 



m(0,h)-m(T, h) 



6 dh \T* 



(17) 

We can calculate the coefficients c x (h) and c m (h) ana- 
lytically in the non-interacting case, 



*4°>(A) 



n 4 ll-3(-Kh/4T*y 



48 \l + (irh/4T*) 2 



(18) 



c<$(h) = 



(h/4T*) 



i + (Trh/4T*y 



(19) 



where T* = ttA/4. 

We can also deduce the values of c x {h) and c m (h) in the 
Kondo regime from the Bethe ansatz results for Xd(0, H). 
As the Wilson or 'x/Y ratio R has the value 2 in this 
regime, independent H as the charge fluctuations can 
be neglected, the value of jd(H) can be deduced and is 
proportional to Xd(0, H). On substituting the value for 
7d(-ff) into equation Ijl6|) c x(^) an d c m (h) can then be 
calculated. The asymptotic values as h — > are c x (0) = 
y/3ir 3 /8 + 0(h 2 ) and c m (h) = hV5TT 3 /16T K + 0(h 3 ). 

In figure 9 we plot the results for — p' d (h) / pd(0) , which 
is proportional to c m (h), for U/ttA = 3.0,0.5,0.0 in the 
range < h/T* < 2.5. It can be seen that all three 
curves have a maximum which implies that c x (h) has a 
zero, and therefore changes sign from positive to negative 
in this range. In the non-interacting case it can be seen 
from the result in equation 118|) that the change of sign 
occurs when h/T* — i/ny/3, and from figure El that it 
occurs for h significantly smaller than T"k in the Kondo 
regime. 



V. LOW TEMPERATURE TRANSPORT IN AN 
ARBITRARY MAGNETIC FIELD 

To extend the calculations to the low energy dynamics 
we use the renormalizcd perturbation theory^, which is 
described briefly in the Appendix. In this form of pertur- 
bation theory we work with the renormalized parameters, 
id, A and U instead of the original bare parameters £d, 
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A and U. The free propagators are those of the non- 
interacting quasiparticles and the expansion is in powers 
of U. As the parameters are already fully renormalized 
counter terms have to be introduced to cancel any further 
renormalization. This expansion is completely specified 
by the three parameters £<i> A and U, and is not re- 
stricted to the low energy regime, but valid for all energy 
scales. In practice, however, the calculations are easier to 
carry out in the low energy regime, where asymptotically 
exact results can be obtained in this regime by working 
only to second order in U. 

Here we exploit the fact that we have these renormal- 
ized parameters as a function of arbitrary magnetic field 
for the symmetric model to calculate the low energy dy- 
namics in the presence of a magnetic field. The quasi- 
particle retarded Green function for the impurity level 
Gd(u>, T, h) takes the form 

1 

G d , CT (w,T, h) 



lo + crid(h) + iA(h) - T, h) 



(20) 

where T, h) is the renormalized self-energy. 

To calculate the leading order T 2 term in the trans- 
port coefficients we need T, a (u>,T) both to order to 2 and 
to order T 2 . We calculate this from the renormalized 
perturbation expansion taken to order U 2 (h). This takes 
full account of the quasiparticle scattering and gives the 
exact result of Yamada 16 for h = 0. The corrections to 
order to 2 can be deduced from the second derivative of 
the self-energy with respect to lo evaluated at to = and 
T = 0. The result for the renormalized self-energy to this 
order is given by 

t a (u, 0, h) = -&(h)co 2 [i - (2 + a u (h))a£(h)] , (21) 

where e(h) = id(h)/ A(h), and a(h) is given by 



a(h) = -p d (0,h)(R(h)-l) 2 , (22) 



and a u (h) by 



a u (h) = -~ 



2I(h)A(h) 



(23) 



where £(h) = Trpd(0, h)A(h)e(h) and 1(h) is the integral 

, did" dtu' 



1(h) 



OO />OC 



G°(o;")G°(^" + ^)[G ( ?(^)] 3 



-oo^~oo 27T 2ir 

(24) 

G^(lu) is the propagator in the diagrammatic expansion 
for T = which is given by 

[G° MP 1 = u + aid(h) + sgn(uj)iA(h). (25) 

The corresponding result for the renormalized self-energy 
to order T 2 can be derived using the Sommerfeld expan- 
sion. The calculation can be performed by using for each 



internal propagator G^(w) in the T = diagrammatic 
expansion an additional correction term, 



(nTf 



S'(u})A{h) 



(u + aed(h)) 2 +A 2 (h) 



(26) 



The result for the renormalized self-energy to order T 2 
for lo = is 

± a (T, 0, h) = -&(h)(irT) 2 [i + ae(h)(l + a T (h))] , (27) 
where the parameter ax(h) is given by 



ar(h) 



1 



6t(h)e(h) 



l-£(/i)tan _1 (e(/i)) 4+^ 



1 



t(h)e(h) 
(28) 

We can now apply these results to the calculation of 
transport coefficients. 



A. Application to magnetic impurities 

The contribution to the conductivity from the scatter- 
ing of isolated impurities described by an AIM, a(T, h), 
given byi^ 



er(T, h) = a ^2 J 



1 



— oo Pd,a(u,T, h) 



duj, (29) 



where pd(u,T,h) = A(h)p~d(u,T, h)/A, and pd(uj,T, h) 
is the spectral density of the quasiparticle Green function 
Gd(co, T, h). The Sommerfeld expansion gives for to 
second order in T, 



a(T,h) = cr { ^2 



Pd,cr(0, T, h) 

PdJOAh) 2 _ p^(0,0,fe) A1 (ttT) 2 ■ 
p d . a (0,0,h) 2 p d , a (0,0,h), 



(30) 



where the prime refers to a derivative with respect to 
lo. Note that the first term still contains a temperature 
dependence via pd- 

On using the renormalized self-energy to calculate the 
quasiparticle spectral density p(u>, T, h), and substituting 
in equation (JJOJ, the final result for a(T, h) to order T 2 



,17 



a(h,T) = a(h,0){l + a 2 (h) (J^tJ + 0(T 4 )}, (31) 
where a(h,0) — 2ao / cos 2 (irm(h)) and a 2 (h) is given by 



<J2(h) 



cos 2 7rm(/i) 



l + C(h)(R(h) -1) 



(32) 



The coefficient C(h) is 

C(h) = 2cos 2 (7rm(/i))-sin 2 (nm(h)) [1 - M T (h) + a^(h)} 

(33) 
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FIG. 10: Field dependent coefficient a% (h) from 
(|32|l for the second order temperature expansion of 
the conductivity. Weak coupling (U/ttA = 0.5) 
up to strong coupling (U/irA = 4) is considered. 



In figure ^3 wc show the second order coefficient 
<72 (fo) plotted over log(/i/T*) for a range of parameters 
(U/ttA = 0.5 — 4). For zero field the conductivity rises 
with temperature as is well known 16 . When h is increased 
02 (/i) decreases and tends to zero for very high fields, so 
that the low temperature conductivity becomes temper- 
ature independent. The impurity level is then shifted out 
of the range of the thermally excited states in the conduc- 
tion band so that there is negligible impurity scattering. 
We note for the strong coupling cases where there is a lo- 
cal moment (U/irA — 2,4), the coefficient 172(h) changes 
sign for a certain critical field h Cl with h c ~ 0.5T*. This 
does not arise, as one could expect from the change of 
the impurity density of states at the Fermi level with the 
field [second term in J5DJ], but from the temperature de- 
pendence of the self-energy, which contributes in the first 
term in (|30|) . To see this only consider the imaginary part 
contribution to the self-energy, 123 which leads to a T 2 
term 

(ttT) 2 /2(R(h) - 1) 2 [cos 2 (ttto(») - sin 2 (7rm(/i))]. (34) 

This expression clearly changes sign when m(h) exceeds 
1/4 or id(h) > A(h). In contrast the contribution from 
the second term in (|30[1 with only the imaginary part of 



the self-energy taken into account gives 

(ttT) 2 /3(1 + (R(h) - l) 2 [cos 2 (7rm(/i)) - sin 2 (7rm(/i))]), 

(35) 

which cannot change sign. Including the real part con- 
tribution does not alter this behaviour substantially. 

Physically, it is the spin flip scattering of the local mo- 
ment that causes the resistance to rise as the temper- 
ature is lowered, leading to a resistance minimum and 
the Kondo effect. Perturbation theory shows that spin- 
flip scattering gives a diverging amplitude for T ~ T*. 
For a mainly polarized impurity spin these processes are, 
however, strongly suppressed, and the change in sign of 
the temperature dependence might be attributed to the 
thermal spin disorder scattering. To our knowledge, this 
effect has not been observed but for magnetic impurities 
systems with a very low Kondo temperature it might be 
feasible to put the result to an experimental test. 



B. Application to quantum dots 

The Kondo effect in a magnetic field has been ob- 
served experimentally in mesoscopic systems, for exam- 
ple, quantum dots in heterostructuresi^. Such systems 
can be quantitatively described by the Anderson impu- 
rity model 19-20 . The symmetric Anderson model cor- 
responds to the situation, where the gate voltage V g is 
tuned to the middle of a Coulomb valley with odd num- 
ber of electrons on the dot. 

Hershfield et aliSi and Meir and Wingreen2£ derived 
an expression for the current through a quantum dot by 
a non-equilibrium calculation, which for the case of sym- 
metric coupling to the leads takes the form 

1 = ^E/^(- ImG r e »)x 

u 

x[n F (uj - Hl) - n F (u - fj, R )], (36) 

where np is the Fermi function, Gj° ncq is the retarded 
non-equilibrium Green function from the Keldysh for- 
malism on the dot, /i^ = \i + eV/2, /i^ = /i — eV/2 are 
the chemical potentials in left, right lead, respectively, 
and V — Vds the source drain voltage^ Go = e 2 /7r?i 
with Planck's constant Ti. 

In the limit of linear response the equilibrium value 
of the one-electron Green function can be used to evalu- 
ate ()36f) . and the resulting expression for the differential 
conductance G — dl/dV through a quantum dot is 

(37) 

In the low temperature regime we can again apply the 
Sommerfeld expansion to obtain the leading order finite 
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temperature corrections to order T 2 , 

G(T,h) = ^[£V Wi<y (0,T,fc) + 



5>/C(o,o,/0 



6 



(38) 



Here, in contrast to the earlier case the second term 
changes sign for a critical field h c when p ( \_ IJ {uj 1 0, h) 
changes from a maximum to a local minimum at u) = 
0. For free quasiparticles this happens when id(h c ) — 
A(/i c )/V3, as discussed in detail in reference 2 ^. The tem- 
perature dependent part of the first term in equation H38|) 
including only the imaginary part from the self-energy 
123 is 



-^(R(h)-l)' 
2A J ' 



, sin 2 (7rm(/i)) — cos 2 (7rm(ft.)) 
(1 + tan 2 (7rm(/i))) 2 ' 



(39) 



and this also changes sign for m(h) > 1/4. The total 
result is 



G(T, h) = G(0, h) 1 



with 



and 



G(0,h) = G cos 2 (7rm(/i)), 



(40) 



(41) 



G 2 (h) 



cos 2 (7rm(/i)) 



— sm' 



{cos 2 (vrm(/i)) [l + 2(R(h) - l) 2 ] 
(irm(h)p + (R(h) - 1) 2 (1 + 2a u (h) - 6a T (h))]\. 



In figure ^2 the field dependence of G2Q1) is shown. 
Note that we have included a minus sign before the T 2 
term in (|40|1 . so that the similar behaviour in figures 
1101 and ITT1 actually corresponds to opposite temperature 
dependence. This is due to the approximate inverse 
relation between the two systems, if the hybridization 
Vk = for an impurity, there is no scattering and hence 
infinite conductivity, whereas if Vk = for the quantum 
dot there is no current and hence infinite resistivity. 




FIG. 11: Field dependent coefficient G2(h) for 
the second order temperature expansion of the con- 
ductance. Weak coupling (U/nA — 0.5) up 
to strong coupling (U/irA = 4) is considered. 



of Bulla et ali 2 ^, in which the self-energy is deduced 
from the calculation of higher order Green functions, as 
this gives more accurate results. In figure 1121 we give 
results for the spin up part of the d-site spectral density 
/Od,|(w) = — i±mGd,t(w + ) for a strong coupling situation 
(U/ttA = 4) for various values of the magnetic field h. 
The shift of the spin-up Kondo resonance from the Fermi 
level with increase of magnetic field, which is almost 
imperceptible on the scale shown, is accompanied by 
large shifts of the spectral weight on the higher energy 
scales as the impurity is magnetically polarized. 



VI. BEYOND THE LOW ENERGY REGIME 

We can use the extension of the NRG method to 
the calculation of dynamic response functions 2 ^* 2 ^ to 
look at the behaviour of the model in an arbitrary 
magnetic field on higher energy scales. In doing so 
it is important to use the density matrix (DM-NRG) 
method introduced by HofstetteflSI as the standard NRG 
approach gives results which considerably underestimate 
the shift of the high energy spectral weight with the 
variation of magnetic field. We also use the approach 



In order to test the improved estimates of this shift from 
the DM-NRG, we calculate (n^ ) by integrating the spec- 
tral density up to the Fermi level to deduce the value 
of the induced magnetization m. This result can then 
be compared with that calculated earlier from the renor- 
malized parameters using equation Jj|J), and also with the 
values calculated directly from the ground state occu- 
pation numbers at the impurity site, as deduced from 
the matrix elements in the NRG routine. The results 
are shown in figure 1131 for the three parameter regimes 
U/ttA — 0.5, 1,4 over ln(/i/7rA). The magnetization de- 
parts from zero already for very low h/irA in the strong 
coupling case due to the induced lower energy scale A, 
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FIG. 12: Strong coupling (U/irA = 4) spectral density 
of the d-site Green function pdj(w) for various mag- 
netic fields h. The energy scale is given by AT* = irA. 



whereas for the intermediate and weak coupling cases 
higher fields are required to show similar behaviour. For 
large h all three merge and approach a saturation value 
77i s = i, corresponding to complete polarization of the 
electron on the impurity site. The agreement of the re- 
sults is excellent in the three parameter regimes for all h, 
and confirms that the spectral weight in a broken symme- 
try situation can be computed correctly using the DM- 
NRG. 
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FIG. 13: Comparison of d-site magnetization m(h) 
for the full range of h for U/irA — 0.5, weak cou- 
pling (w), U/ttA = 1, intermediate coupling (m) and 
U /ttA — 4, strong coupling (s) computed from the 
renormalized parameters (RP), the occupancy (occ) 
and the spectral weight as described in the text. 

In figure we focus on the shift of the quasiparticle 
(Kondo) resonance in the results for the strong coupling 
case (U/ttA = 4) as in figure IT21 This shift of the reso- 
nance from the Fermi level (ui = 0) with increasing mag- 



netic field values is clearly seen on the higher resolution 
energy scale used for this plot. As the peak shifts, its 
height decreases and the resonance becomes broader. For 
even larger fields than shown here the peak merges with 
the lower atomic limit peak seen in figure El Note that 
the peak form is asymmetric with logarithmic tails, sim- 
ilar to the results of Rosch et obtained using the 
perturbative RG for the Kondo model. However, some 
of the asymmetry in the our results must be attributed 
to the logarithmic broadening scheme. 
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FIG. 14: Quasi particle peak for the spectral density 
of the d-site Green function p^(u>) for various values 
of h for the same parameter set as used for figure ITS! 

If —e p (h) is the position of the quasiparticle peak in the 
spectral density for a spin up electron, then the corre- 
sponding value for non- interacting electrons (U = 0) is 
half the Zeeman splitting, Az, where Az — 2h. An ex- 
act expression for e p {h)/h in the limit h — > has been 
derived by Logan and Dickens^, 



lira 

h->0 



e p (h) 



R 



1 + bAz 2 ' 



(42) 



where b is the curvature of the imaginary part of the self- 
energy at u> = 0. The value of b can be calculated from 
the renormalized perturbation expansion 4 ! and the result 
(I42|l written as 



lim 

h— >o 



£ P (h) 



R 



1 + (R- l) 2 /2 



(43) 



This ratio, therefore, varies from 1 in the non-interacting 
case (R = 1) to 4/3 in the Kondo limit (R = 2). It is 
not straight forward to obtain a precise estimate of b or 
the value of e p (h) from the DM-NRG spectra as they are 
sensitive to parameters of the logarithmic scale Gaussian 
broadening which is used to obtain a continuous spec- 
trum on all energy scales from the discrete results. How- 
ever, if the broadening is modified to Lorentzian peaks 
with constant width for the very low energy scales the 
asymptotic results can be confirmed. 
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We have estimated the ratio e p (h)/h from the NRG 
spectra for higher magnetic field values and find that the 
it increases monotonically with ft and exceeds the value 
of 2 before the peak merges at high field values into the 
atomic-like peaks. There have been other estimates of 
the /i-dependence of this ratioSi^i*^, but these differ 
markedly according to the method of calculation. On 
the basis of a Bethe ansatz calculation of the spinon 
spectrum for the Kondo model Moore and WenSi find 
that £ p (h)/h < 2 in all cases and conjecture that the 
value of 2 is the high field asymptotic limit. It is possible 
that this is a feature of the localized model, when charge 
fluctuations are completely suppressed. There is some 
evidence in support of this in our results in that, as 
we suppress the charge fluctuations on increasing the 
value of U through the values U/wA = 2,3,4, the ratio 
e p {h)/h increases less rapidly with increase of ft. The 
ratio only begin to exceed the value of 2 roughly at the 
point when charge fluctuations set in and i?(ft) begins 
to differ significantly from the value of R(h) for the 
localized model, R(h) = 2. Costi^ has also done NRG 
calculations for a localized model and finds a ratio close 
to but always less than 2. Using the local moment 
approximation Dickens and Logan^S have also estimated 
the ratio e p {h)/h and find an even more marked increase 
in the ratio with increase of ft to values such that 
e p (h)/h > 2. 

We have also studied the dependence of the ratio 
e p (h)/h on the Kondo temperature Tr. In figure 1151 we 
plot the results starting from a minimal Tk,o correspond- 
ing to U/irA — 6 for various fields ft. One finds that 
with decreasing Kondo temperature, which translates to 
increasing renormalization of the quasiparticles, the peak 
splitting and thus the sensitivity towards the exposure to 
a magnetic field is enhanced. The results can be shown 
to be fitted well by a logarithmic law, 



e p (h, Tk) _ e p (h, Tk,o) 
ft ~ ft 



a In 



Tk, o 



(44) 



For all the curves one finds a = 0.26 ± 0.01 independent 
of ft. 

For the total spectral density, which we denote by pd(uj), 
we need to include the contribution from the down spin, 
which has its peak at £ P (h), pd = Pd.T+Pd,!- The splitting 
between the up and down peaks in the total spectrum is 
then 2s p (h)f c (h), where f c (h) is a correction factor due 
to the overlap of the resonances^ In general it has to be 
determined numerically, but for free quasiparticles (see 
equation l|20|) without it is given by 



f c {hf = i - i - 



A(ft)/e d (ft) 



1/2N 



(45) 



For higher fields one has f c (fjb,A) ~ 1 — A 4 /8(?76) 4 . In 
figure lT^l we give an example for the total spectral density 
for the earlier used parameters U /ttA — 4 and a range of 
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FIG. 15: The dependence of the ratio e p (h)/h 
on Tk for various ft. The energy scale 

T* is given by the corresponding values for 
U/nA = 4, which satisfies T*/T K , ~ H- 



magnetic fields. Clearly, the peak splits above a critical 
field, h > 0.5T*, which is in agreement with results of 
Costi for the Kondo model^. 
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FIG. 16: Total spectral density of the d-site pd(^) 
for various fields ft. One can see that the peak 
splitting becomes visible only for fields h > 0.5T*. 



Comparison with experiments on quantum dots 

In section V.B we have calculated the lowest order 
temperature dependence of the conductance through a 
quantum dot in a magnetic field, (|40|l . This temperature 
dependence has been observed in the zero magnetic field 
case— i. We noted earlier that there is a sign change in 
this leading temperature dependence at a values of the 
magnetic field < ft < T*. The sign change in the sec- 
ond term in equation (|37|) occurs when pd changes from a 
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local maximum to a minimum as can be seen in figure [TBI 
A qualitative explanation of this sign change is that the 
local spectral density on at the Fermi level is suppressed 
with increasing magnetic field. At higher fields when the 
spectral density develops two peaks then there are more 
thermally excited states which can contribute to an in- 
crease of the conductance. This temperature dependence 
might be experimentally observable, since estimates of 
the Kondo temperature are of the order 300mK corre- 
sponding to magnetic fields in the experimental ranged 
A difficulty might be that the overall response is reduced 
by the cos 2 (7rm(/i)) factor in equation (|4lip. 

Some of the quantum dot experiments 3 ^* in the pres- 
ence of a magnetic field have been performed in non- 
equilibrium situations with a finite source-drain voltage 
V. Two peaks can be observed in the differential conduc- 
tance as a function of the voltage V for fields strengths 
larger than a critical value. There have been several 
interpretation s 30 ! 31 of these results based on the approx- 
imation of using the equilibrium Green function to eval- 
uate 1(36(1 . With this approximation at T = we get 
an expression for the differential conductance G(V) as a 
function of the voltage V, 

G(V) = § = ^ Pd (eV/2). (46) 

In this approximation G(V) is directly proportional to 
the total equilibrium spectral density evaluated at u> = 
eV/2, which is shown in figure [TBI To test whether the ex- 
perimental results can be explained on the basis of equa- 
tion 1)46(1 we have extracted the magnitude of the peak 
splitting for U/nA = 2,4 and compared it with experi- 
mental results*^, which are displayed in figure ITTI 

We can see in figure El that, whilst there is an agree- 
ment in the range /i/Tk — 1, in general there does not 
appear to be a satisfactory explanation of the experimen- 
tal results based on approximating the non-equilibrium 
Green function by the equilibrium one. Hewson, Bauer 
and OguriSa have calculated the non-equilibrium Green 
function for a small but finite voltage V, and find that 
a finite voltage reduces the peak position substantially. 
We conclude that an accurate agreement of experiment 
and theory rests therefore on an accurate description of 
the steady state situation out of equilibrium. We are 
currently working on an extension of renormalized per- 
turbation theory in the non-equilibrium formalism, which 
takes this situation into account appropriately. 



VII. SUMMARY 

We have shown that in an arbitrary magnetic field H 
it is possible to describe the quasiparticles of a Fermi 
liquid regime by field dependent parameters. We have 
used the particle-hole symmetric Anderson model as an 
illustration because of its wide range of applications and 
because it serves as a local model of strong correlation 
physics. For this model the three relevant field dependent 
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FIG. 17: Peak position V p of the theoretical results 
from 1)46(1 obtained from the NRG spectra in compari- 
son with experimental result Ak (referred to as half the 
Kondo splitting) 33 . The experimental Kondo tempera- 
ture is inferred from the critical field for peak splitting 
B<° xp) « 2T and the result h c « 0.5842T K for strong 
coupling^ A dashed line corresponding to twice the Zee- 
man energy has been included for comparison. Note that, 
although all curves lie above 2h, due to u> = V/2 in the 
argument of pd in equation 1(461) the actual peak splitting 
does not exceed twice the Zeeman splitting in all cases. 35 



parameters are the renormalized impurity levels with spin 
er, £(i jCT (i?), the quasiparticle resonance width, A(H) and 
the local interaction U{H) . We have shown how these can 
be deduced from NRG calculations of the low lying ex- 
citations. Once the renormalized parameters are known, 
the impurity spin and charge susceptibility, the specific 
heat coefficient and the induced impurity magnetization 
at T = for arbitrary magnetic field can be expressed ex- 
actly by substituting into the relevant formulae derived 
from a renormalized perturbation theory. The leading 
temperature dependent corrections to the susceptibility 
and magnetization can be obtained in a similar way. We 
have also extended the renormalized perturbation expan- 
sion in order to calculate the leading temperature depen- 
dence for the finite conductivity due to scattering from 
an impurity in a metallic host, and for the conductance 
through a quantum dot. 

By choosing the bare parameters in the absence of the 
field to correspond to the strongly correlated or Kondo 
limit, the de-renormalization of the quasiparticles can be 
followed as the magnetic field strength is increased from 
zero. For extremely large magnetic fields the parame- 
ters revert to their bare values. This approach gives an 
overview of the low energy behaviour of the model as a 
function of the applied magnetic field strength. 
A number of physical properties are found to change 
qualitatively in the strongly correlated case for magnetic 
field strengths in the range < g^H < 7k, where 7k 
is the Kondo temperature. This should be a physically 
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accessible magnetic field range for many systems. The 
T 2 coefficient of the magnetic susceptibility, the conduc- 
tivity from a magnetic impurity in the strong correlation 
regime, and the conductance through a quantum dot all 
change sign in this magnetic field range. 

The approach developed here is a general one and is 
equally applicable to asymmetric impurity models and to 
lattice models. For lattice models, for which dynamical 
mean field theory is applicable, similar NRG methods to 
those employed here can be used^, and our calculations 
are currently being extended to models of heavy fermions. 
We note that the approach is not restricted to the NRG 
method, the relevant renormalized parameters could also 
be estimated using other theoretical techniques, varia- 
tional methods for example. 
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VIII. APPENDIX: THE RENORMALIZED 
PERTURBATION APPROACH 

We give a brief synopsis of the renormalized pertur- 
bation approach. This approach is best developed in 
the field theoretical Lagrangian formalism where the 
renormalization of the field is a more natural concept. 
The Lagrangian corresponding to the bare Anderson 
£AM(£d, A, U) can be rewritten in the form 

£AM(£d,er, A, U) = ^AM^d.cr, A, U) + £ ct (Ai, A 2 , A 3 ), 

(47) 

where the renormalized parameters, id. a- and A CT , are de- 
fined in terms of the self-energy S CT (lu) of the one-electron 
Green function for the impurity state, 

G^M = * v ( r , (48) 



and are given by 

ed ) <7 = z CT (£ d) <r + S CT (0)), A CT = z (T A, (49) 

where z a is given by z a = 1/(1 — S^(0)). The renormal- 
ized or quasiparticle interaction U, is defined in terms of 
the local irreducible 4-vertex r-|-j_ (a;i , U)2, ui^, U4) at zero 
frequency, 

L> = z T z i r n (0, 0,0,0). (50) 

The remaining part of the Lagrangian £ ct (Ai, A2, A3) con- 
tains three counter terms with coefficients Ai, A2 and A3. 
A perturbation expansion in powers of the quasiparti- 
cle interaction U can then be carried out provided the 
counter terms are also taken into account. The counter 
terms Ai, A2 and A3 are determined to each order in U 
by the renormalization conditions*, 

£<r(0)=0, E' CT (0) = 0, f u (0,0,0,0) = l> (51) 

where and T^i(u>i, u>2, W3, W4) are the renormalized 

self-energy and irreducible 4-vertex, respectively. 

The renormalized parameters, as defined in equations 
(|49|l and 1)50(1 . can be calculated in low order perturba- 
tion theory*, but these calculations cannot be extended 
to the more physically interesting strong coupling regime 
U/ttA > 2. These parameters, however, can be identified 
with those in equation @ , and so they can alternatively 
be deduced from the levels in the numerical renormal- 
ization group calculations, which can be carried out for 
all values of U. The expressions given earlier for the im- 
purity occupation number at T = 0, the specific heat 
coefficient and spin and charge susceptibilities, ©, © 
and JBJ, are exact and correspond to first order calcula- 
tions in a renormalized perturbation expansion* in pow- 
ers of U. Calculations carried out to second order in U 
give exact results for the leading low order temperature 
dependence to the impurity resistivity* and the leading 
non-linear term in the differential conductivity**^*. 
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